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Abstract 



Q^ We look at the response of a nonlinearly coupled scalar field in an asymptotically AdS 



black brane geometry and find a behaviour very similar to that of known dissipative non- 
linear systems like the chaotic pendulum. Transition to chaos proceeds through a series 
of period-doubling bifurcations. The presence of dissipation, crucial to this behaviour, 
J> arises naturally in a black hole background from the ingoing conditions imposed at the 

On horizon. AdS /CYT translates our solution to a chaotic response of O, the operator dual 

^ to the scalar field. Our setup can also be used to study quench-like behaviour in strongly 

\Q coupled nonlinear systems. 
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1 Introduction 

Over the last few years AdS /CFT [1-4] has been successfully applied to various problems 
inspired by condensed matter systems [5]. One kind of system which is of much interest, is 
one which is driven or quenched externally by a forcing [6]. In the field theoretic context it 
translates to creating a disturbance by turning on a time dependent source for some operator. 
In general, in a strongly coupled field theory, it is difficult to study the driven problems directly 
by employing analytic or even numerical methods. Holography gives an alternative avenue to 
this strongly coupled physics. Holography maps strongly coupled field theoretic problems 
to gravitational problems involving Einstein's equations and additional equations governing 
the dynamics of dual bulk fields. The resulting differential equations, fewer in number, can 
be studied numerically or analytically. Various authors have discussed the issue of driven 
systems in the context of holography [7-13]. For example, the problem of slow quench has 
been discussed in the simplest of setups in [10]. In this paper, we discuss the problem of 
nonlinear dynamics in the same spirit in context of holography. 

We start with a simple model of a scalar field in an AdS black-hole geometry [14]. We drive 
the system by introducing a time-dependent source for the scalar operator on the boundary. 
Without any nonlinearity, the problem boils down to the well known calculation of scalar two 
point function in AdS black hole [15]. Here we would like to understand how nonlinearity 
changes this simplistic behaviour. There can be two different ways to approach this problem; 
one is to consider the gravitational back reaction [9] of the scalar field and the other is to 
turn on the scalar self interaction. Here we choose the second approach and consider the 
scalar field as a probe field with negligible gravitational back reaction. One can think of the 
scalar as coming from some brane field. The question we ask here is what kind of steady 
state behaviour, if any at all, do we observe. Our system can be thought of as a driven 
diffusive system like a forced damped pendulum. In our case, absorption at the black hole 
horizon introduces a diffusive element. Not unexpectedly, at small values of the nonlinearity 
parameter the effect of nonlinearity is negligible and we get a steady state solution. However, 
as we increase nonlinearity gradually we see period doubling and novel solutions in bulk with 
response at frequencies lower of that of the driving frequency. As we increase nonlinearity 
more signatures of nonlinear dynamics are seen as chaotic motion sets in^. 

Chaos, a behaviour seen in many nonlinear systems, lacks a precise definition but can be 
roughly be defined as an exponential sensitivity of the dynamics to its initial conditions [21]. 
Chaos is more frequently characterized by the qualitative behaviour of the system in the chaotic 
regime (which on the other hand is quite independent of the initial conditions), or by the route 
the system takes to chaos as a control nonlinearity parameter is varied. In various complex 
systems, many different routes to chaos can be seen - period-doubling bifurcations, quasiperi- 
odicity, intermittency etc. Quasiperiodicity and intermittency are characteristic of systems 
with multiple natural frequencies and systems with higher degrees of freedom respectively. We 
will restrict our attention to the period-doubling route of transition to chaos. 

The rest of this paper is organized as follows. In Section 2, we review some of the features 
seen in simple nonlinear dissipative systems by taking an example of the forced damped pen- 



^For other examples of chaotic dynamics in the context of string theory see [16-20]. 



dulum. In Section 3, we set our problem up and obtain the equations of motion. We solve the 
equations of motion numerically and in Section 4 we display the results for the static case and 
with time-dependent boundary conditions relevant to both chaotic behaviour and quenches. 

2 Chaos in the driven damped pendulum 

In this section, we review some of the features of the response seen in nonlinear dissipative 
systems. For an extensive discussion we refer the reader to numerous standard textbooks 
available [21, 22]. The reader familiar with period doubling bifurcations and related aspects of 
dynamical systems may proceed directly to Section 3. A classic example of the above behaviour 
is seen in the forced damped pendulum described by the equation of motion, 

e + -fO + n^ sine = F{t). (i) 

Here Q is the natural frequency of the pendulum (for a simple pendulum of length / oscillating 
due to gravity, Q, = \Jg/l) and 7 is the damping coefficient. We are not restricted to small 
amplitudes and the angle 9 can be large - a nonlinearity will arise due to the presence of 
the sin 6' term. The driving is typically given by sinusoidal forcing with amplitude A and 
frequency Vt: 

F{t) = Asinnt. (2) 

For short times, there is a usually transient behaviour, which will not be of concern to us here. 
We will look at the behaviour at late times after the transience has died down. This behaviour 
is qualitatively different for different driving parameters [A^Vt). For small amplitudes A, the 
late time behaviour is insensitive to the precise initial conditions we begin with. All trajectories 
in the phase space converge to a periodic attractor or a limit cycle. For very small amplitudes, 
the response frequency is at the forcing frequency Q and is roughly sinusoidal. The phase 
portrait is a closed curve, roughly an ellipse - we will call this limit cycle a 1-cycle. As the 
forcing amplitude is increased, the phase space ellipse distorts. However at a particular value 
of the amplitude the ellipse abruptly splits into two - the alternate crests and troughs begin 
to have slightly different heights. The response is thus no longer at frequency Q now, but 
Q/2. This transition is known as period- doubling. As mentioned before, the phase portrait 
for the motion is now not one but two intersecting loops - this limit cycle is a 2-cycle. As 
the forcing amplitude is increased further, we get subsequent period-doubling more and more 
rapidly to 4-cycles [Fig. 1(a)], 8-cycles and so on. We eventually reach a regime where there is 
no periodicity at all and there is no closed limit cycle. In this regime the response is chaotic - 
the aperiodicity is accompanied by the fact the system is exponentially sensitive to its initial 
condition. An efficient way to represent the period-doubling transitions is to plot what is known 
as the bifurcation diagram [Fig. 1(b)]. For each value of the control parameter (amplitude here), 
one looks at all the values of ^ or ^ attained after intervals of multiples of the forcing period - so 
this is a stroboscopic sampling of the phase-space. For a 1-cycle only one value is attained, for 
a 2-cycle two different values, for a 4-cycle four values etc. and the period-doubling transition 
from one to the other is clearly seen as a splitting or a bifurcation in this plot. 



The period-doubling bifurcation route of transition to chaos is generic to nonhnear systems 
with low number of degrees of freedom including even very simple examples like the tent 
or logistic map [21, 22]. The existence of period-doubling bifurcations can for all practical 
purposes be treated as a positive signature of chaotic behaviour in the system. Some other 
positive indicators of chaos are the presence of 3 and other odd cycles past the chaotic regime, 
which we will also see in our system below. 

With this very quick introduction to chaos in non-Hamiltonian dynamical systems, we move 
on to describe our system, a nonlinearly coupled scalar field in a background AdS black hole 
geometry. 
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(a) Phase portrait showing 4-cycle. 



1.065 1.070 1.075 1.080 1.0 

(b) Bifurcation diagram. 



Figure 1: Nonlinear response of the pendulum (with ^2 = 1, 7 = 0.5). Fig. 1(a) shows the phase 
portrait for Q = 2/3, A = 1.081. We see four lobes corresponding to a 4-cycle. Fig.l(b) shows 
the bifurcation diagram, typical to this route of transition to chaos. 



3 Setup 



We consider a black brane geometry in asymptotically AdSd+i described by the metric: 

W) 

where the warp factor 



ds^ = -f{r) dt^ + ^ + r^ dx^ 



(d-i) > 



/(r)=rMl-J 



(3) 



(4) 



In these coordinates, the AdS boundary is located at r = 00 and the black brane horizon at 
r = tq. In order to preserve regularity at the horizon, we work in ingoing Eddington-Finkelstein 
coordinates, 

dr 



v = t + 



fir) 



(5) 



in which the metric becomes 

ds"^ = -f{r) dv^ + 2dvdr + r^ ^x^d-i) • (6) 

We add a minimally coupled scalar field $ with a nonlinear potential term, given by the 
Lagrangian : 

S = -J rf"+'xv^ I ^(7»^V„$V,$ + ^m2$2 + iA$^ I . (7) 

We work in the probe limif* where the field $ does not backreact on the geometry [14]. The 
$^ interaction leads to a nonlinear term in the equation of motion, which is given by 

V^$-m2$- A<l>^ = 0. (8) 

Here 



25,5,$ + f{r)d^,^ + ^-^5„$ + li^^UM. + f'(r)\dr<^ . 

r I r J 



It will be useful to make a coordinate transformation 



In d = 4, which will be of relevance to us. 



1 f n . -1 f ^ \ 1 f^~fh 



d-l 

r ' 2 



(9) 



dv 
P= I jT-^ , such that - f{r)dr = dp , (10) 



pW^-|.-2ta„-'^-J-log^^J|. ,11) 

The boundary is now mapped to p = and the horizon to p = oo. V^$ simplifies to 

V2$ = -^ l-2d,dp^ + d^^ + ^'^'^^^^''^ (5,$ + 5,$) I . (12) 

After a field redefinition [25], 

$ = r-'^(/), (13) 

V^$ simplifies even further to 



V^'l' = j^ {-2d,dp<P + dl<P - r(r)0} , (14) 



*It should be noted that in the non-probe case there is always a horizon formation [23, 24]. 



with 

^^, ^_ {d-l){d-3)p{r) , id-l)fir)f'{r) 

Then the full equation of motion takes the form, 

2dydp(j) - dlcj) + [V{r) + m^fir)] (j) + Xr-^'^-^^ f{r)(j)^ = . (16) 

Close to the boundary, where /(r) ^ r^ ^ p^^, the two solutions scale as 

(pip ^ 0) ~ p^^i , with A± = ^-— ^ . (17) 

We will choose to work with w? = (1 — (i^)/4 which is the conformal mass (it is still above the 
Breitenlohner-Freedman [26, 27] bound — (i^/4). For this value of ?7i^, A-t = 1,0 and 

0(p ^ 0) ~ 00 - P01 + • • • (18) 

Near the boundary, the bulk field $ scales as: 

^{p,x^) = p'-^ {p'^-M^n-p^^M^n} ■ (19) 

Using the standard AdS /CFT prescription, one interprets 0i as the expectation value (O) of 
the operator dual to the field $ sourced by its boundary value J = 0o- Since in our case 
A+ — A_ = 1, (O) is given simply by 



^^)-^--i 



(20) 

p=0 



We need to augment our equation of motion (16) with appropriate boundary conditions. In- 
going condition at the horizon is ensured by regularity of at p = oo in ingoing coordinates. 
For the condition at the AdS boundary at p = 0, we choose a time-dependent "driving" 

*,.ow=^o(»)={;^^::;^;„ (21) 

Existence of a black hole horizon adds the element of absorption in our setup. On one side 
energy is pumped into the system by time dependent boundary condition and on the other 
side energy is absorbed at the horizon. The effect of nonlinearity is to mix up the modes in 
the bulk of the region [Fig. 2]. 

As we will see, using the different forms of driving, our system can be used to study forced 
damped dynamics as well relaxation in quenches. We will look at the various cases in the next 
section. 
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Figure 2: Schematic diagram showing energy flow from boundary to horizon. In the bulk, the 
energy flow is affected by nonlinearity due to mixing up of modes. 

4 Numerical solution and results 

We solve the equation of motion (16) using the method of lines. We discretize on a uniform 
grid in p. The regularity condition at p = oo is identical to a Neumann condition dp(j)\ ^^ = 
and this implies that (p attains a constant value for large p. We can thus cut off the domain 
of our problem to a large but finite value of p. It serves our purpose to choose pmax = 5.5 and 
A^ = 2500 points on the grid.^ We solve the resulting ODEs, which are N coupled first-order 
differential equations in v, using the SUNDIALS IDA solver [28]. 

For the purpose of numerics, we choose d = A and A = 1. However any non-zero A can be 
scaled to A = 1 using A2 — > 0. The forcing amplitude A thus also serves as a nonlinearity 
parameter. Before going over to the dynamics of the system, we will briefly look at the time- 
independent solution. 



4.1 Static case 

For the time-independent case, the equations are ODEs with mixed boundary conditions and 
the solution can be simply obtained, for example, by shooting. The results are shown in Figure 
3. A typical profile of the bulk field is plotted in Fig.3(a). Fig. 3(b) shows the response of the 
field theory operator expectation {O) = —d(f)/dp as a function of the boundary value 0o- We 
see that for large values of the field (high nonlinearity), the operator scales as {O) ~ (pQ . 

This exponent can also be obtained analytically by scaling arguments near the boundary. 
For small p, the potential and mass terms cancel each other for m"^ = (d^ — 1)/4 and the static 
equation is 



-dl<f>{p) + Xp'~'<P'{p)=0. 
If (psip) is a solution to this equation, then 0(p) = a^^'^^^'^cpsiap) is also a solution. 



(22) 



(23) 



^It is clear from Figs.3(a)-4(a) that (j){p) does indeed go to a constant for the Pmax we choose. 
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(a) Typical profile of bulk field. 



(b) Static response of boundary operator. 



Figure 3: Static solution. Fig.3(a) shows a typical profile of the time-independent bulk field. 
Fig. 3(b) is a logarithmic plot of (O) vs. boundary value 0o- {O) ~ 0o ^^^ large 0o- 

Thus (O) ~ 00 . The near boundary scaling works because for large boundary values 

00 the field drops off faster with p and its behaviour is dominated is by the small p solution. 



4.2 Sinusoidal driving 

We first choose a driving 0o(f) = AsinQv. In the limit A ^ 1, the effect of nonlinearity is 
negligible and we have the relation 



{0){n) = G{Q)A, 



(24) 



where G{Q) is the Green's function in the black hole background. 

As we turn on nonlinearity parameter the response (O) ceases to be sinusoidal and a linear 
function of amplitude. However the response is still periodic. We get novel nonlinear steady 
state solutions. Depending on the initial conditions, there is a transient behaviour, which at 
late times usually settles into a steady limit cycle. The results are shown in Figures 4 and 5. 
For Q = 14.1 and A = 7.0, for example, the response is at the driving frequency - the limit 
cycle is a 1-cycle. Typical profiles of the bulk field when the boundary value crosses zero (for 
the late time response) is shown in Fig.4(a). The time-dependent response of the operator (O) 
is plotted for two parameter values in Fig.4(b). 

As the amplitude A is increased keeping the frequency Q fixed, we see a change in the 
qualitative behaviour of the late-time response. For the same frequency Q = 14.1, if the 
amplitude is increased to A = 8.25, the alternate crests and troughs have different heights 
- we see a 2-cycle. Indications of this are also seen in Fig.4(a), where we see two different 
field profiles for alternate periods. The phase portrait for these parameters [Fig.5(a)] has two 
distinct lobes clearly showing a 2-cycle. The system has thus undergone a period-doubling 
bifurcation between A = 7.0 and A = 8.25. From the bifurcation diagram [Fig.5(b)], as the 
amplitude is increased keeping Q = 14.1, we see a clear bifurcation near A ^ 8.0. For larger 





(a) Late-time bulk field (for v = 2n7r/J7). 



(b) Late-time behaviom' of {0){v). 



Figure 4: Sinusoidal driving. In Fig.4(a) we see the late-time profile of the bulk field sampled 
after every period. Two distinct profiles indicate that the field takes two different configurations 
which repeat after every two periods. In Fig.4(b) we look at the time-dependence of {O). The 
driving (upper blue) is at fi = 14.1. For a smaller amplitude A = 7.0 (middle red), the 
response is at the driving frequency (1-cycle). For a larger amplitude A = 8.25 (lower green), 
the response is at half the driving frequency (2-cycle). 

values of amplitude A ^ 9 we see a transition to complete chaos. As we increase the amplitude 
A further, we end up again in a state which is "less" chaotic and looks quasiperiodic [Fig. 6]. We 
also find approximate odd periods. This not unlike what we see in a simple damped pendulum. 
Our system is a 1 -|- 1 dimensional system with an infinite number of degrees of freedom; it is 
unlikely that our chaos to order transition would be sharp and clean. Also one can question 
whether we ever get to an exact strange attractor like configuration. It seems that the entire 
configuration space is explored once the system becomes chaotic. Further resolution of these 
questions need more intensive numerics. 

It should be noted that the value of A at which chaotic motion sets in varies inversely with 
Q. At small value of ^2 <^ 1, the system shows adiabatic behaviour even for large A. At small 
Q, reverse is true and chaotic motion sets in for smaller A. 



4.3 Tanh driving 

We now choose (poiv) = Ata.nh.(3v. This form of the driving can be used to study quantum 
quenches using holography (see, for example, [10]). Our results are shown in Figure 7. For 
a very slow driving /3 <^ 1, the response of the operator {O) is quasi-static - the expectation 
value at any time is same as the expectation value for the same 0o for the static case. For 
a slow variation (/3 =5), there is a slight lag but (O) essentially follows its quasi-static value 
(C)qs [Fig.7(a)]. For a faster variation {(3 =200), we see an oscillatory behaviour about the 
quasi-static value in the response [Fig. 7(b)]. Our results have some of the features of the results 
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(a) Phase portrait showing a 2-cycle. 
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(b) Bifurcation diagram. 



Figure 5: Transition to chaos. In Fig.5(a), the phase portrait has two lobes indicating a 2-cycle. 
In Fig. 5(b), we see a period doubling bifurcation and transition to chaos. 
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Figure 6: Transition from chaos back to approximate order as the value of the amplitude is 
increased. For A = 8, we see a chaotic behaviour (middle red). For A = 9.75, we see a roughly 
periodic behaviour again. 

of [10], however there are crucial differences due to the presence of dissipation in our system.^ 
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